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ABSTRACT 

We perform Monte Carlo simulations to study E p -Ei so correlation in the context of a multiple subjet model (or 
inhomogeneous jet model) for gamma-ray bursts (GRBs), X-ray-rich GRBs (XRRs), and X-ray flashes (XRFs). 
For a single subjet, we find that E p oc Ei SO 0A for large viewing angles. For the multiple subjet model in which all 
the subjets have the same intrinsic properties, off-axis events show E p oc E lso a with 0.4 < a < 0.5. If the intrinsic 
properties of the subjets are distributed so that on-axis emission of each subjet follows a correlation E p oc L^o > 
we obtain the Amati correlation (E p oc Ei SO x l 2 ) over three orders of magnitude in E p . Although the scatter around 
the Amati correlation is large in the simulation, the results are consistent with the observed properties of GRBs 
with known redshifts and the BASTE GRBs with pseudo redshifts derived from the lag-luminosity correlation. 
We also calculate the event rates, the redshift distributions, and the Tyo duration distributions of GRBs, XRRs, 
and XRFs which can be detected by HETE-2, assuming that the source redshift distribution is in proportion to 
the cosmic star formation rate. It is found that the event rates of three classes are comparable, that the average 
redshift of the XRRs is a little larger than those of the GRBs and the XRFs, and that short XRRs arise when a 
single subjet is viewed off-axis or viewed on-axis with slightly high redshift. 
Subject headings: gamma rays: bursts — gamma rays: theory 



1. INTRODUCTION 

HETE-2 observations have provided strong evidence that 
softer and dimmer gamma-ray bursts (GRBs) smoothly ex- 
tend to X-ray flashes (XRFs) through an intermediate class 
of events called X-ray-rich GRBs (XRRs). For events 
with known redshifts and well observed spectra, the rest- 
frame spectral peak energy E p and the "bolometric" isotropic - 
equivalent 7-ray energy Ei SO have a strong correlation, i.e., 
E p oc Ei SO x / 2 (lAmati et alJ 120021) . This E p -Ej S0 correlation, 
called the Amati correlation, has recently b een exte nded down 
to lower energies characteristic of XRFs (Lamb et al. 2004). 
Since various observed quantities other than the Amati corre- 
lation also distribute co ntinuo usly among GRBs, XRRs, and 
XRFs ( Sakamo to et al.l 12005 ). it is strongly suggested that 
these three classes are related phenomena. 

While m any different models have been pro posed for 
XRFs (see iGranot. Ramirez-Ruiz. & Pernal 120051 and ref- 
erences therein), we have proposed the "off-axis model" 
J Yamazaki. Ioka. & Nakamura! 120021 120031) in which XRFs 
are the usual GRB jets viewed from an off-axis viewing an- 
gle (see also Woods & Loeb 1999). When the jet is ob- 
served off-axis, the emitted photons are out of the beam- 
ing cone and less blueshifted than photons emitted along 
the jet axis, so that the events look like XRFs. It has been 
shown that the viewing angle of the jet is the key param- 
eter to understand the various properties of the GRBs and 
that the luminosity-variability/lag/width correlations might 
be naturally derived in the framework of off-axis models 
Jloka & Nakamural2001l) . 

As fo r the Amati correlation. lYamazaki. Ioka. & Nak amura 
( 2004a) computed E p and Et so using the uniform jet model and 
found that the results are compatible with the observations. 
They also found tha t E„ oc Ei SO l ^ in the smaller Ej m regime. 
Eichl er & Levinsonl ( 120041) investigated the correlation in an 
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annular jet model, and derived that if the viewing angles are 
within the annulus, E p oc Ej SO a with 1/3 < a < 1/2, which is 
compatible with the observations. Compared with our uni- 
form jet model, in the annular jet model the energy is large 
due to the emissions fro m widely d istributed segm ents with 
similar viewing angles. Eichler & Levinson (2004) also an- 
ticipated that multiple discrete emissions could have the same 
effect. 

The off-axis jet model has recently been improved to in- 
clude short GRBs (Yamazaki, Ioka, & Nakamura 2004b) as a 
unified model, where the GRB jet is not uniform but made 
up of multiple subjets or multiple emission patch es. This is 
an ex t reme case of an inh omogeneous jet model ( Na kamural 
2000; Kumar & Piran 2000). The crucial parameter is the 
multiplicity, n s , of the subjets along the line of sight. If 
n s > 2, the burst looks like a long GRB, and if n s = 1 the 
burst looks like a short GRB, while if n s = the burst is an 
off-axis event for all the subjets and looks like an XRF or 
an XRR. We also found that the unified model may explain 
the bimodal distribution of the 7qq d urations of BATSE GRBs 
( iToma. Yamazaki. & Nakamural2005l) . 

In this paper, we examine E p -Ej S0 correlation in the multi- 
ple subjet model to show that the unified model is consistent 
with the observations of E p and Ej S0 . This paper is organized 
as follows. In §|2 we describe our multiple subjet model for 
prompt emissions. First, the E p -Ej SO correlation for a single 
subjet is discussed in § [3] and then we discuss the results of 
Monte Carlo simulations in the multiple subjet model in § 
Section|5]is devoted to discussion. 

2. PROMPT EMISSION MODEL 

Let us suppose that iV tot subjets with opening half-angle 
ub ^ l auncne d from the central engine of GRB ran- 
domly in time and directions and that the whole jet with open- 
ing half-angle A6 tot consists of these subjets. We introduce 
the spherical coordinate system (r,-d,<p) in the central engine 
frame, where the origin is the location of the central engine, 
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and -& = is the axis of the whole jet. The axis of the jth subjet 
(J = 1, • • • ,Mot) is denoted by (i9 u \ ip u> ). If the direction of the 
observer is given by (# bs) fobs), the viewing angle of the jth 
subjet from the line of sight is 

9f = cos -1 [sin $ b s sin $ tJ) cos(<^ b S - f' J) ) + cos i9 Q b s cos $ 0) ] . 

(1) 

For emission model of each subje t, we use the same formu- 
lations and notations as used in Yamazaki, Ioka, & Nakamura 



(2003). Let us use another spherical coordinate system 
(r, 9, <f>) in the central engine frame, where the origin is the lo- 
cation of the central engine, and 9 = is the line of sight. We 
adopt an instantaneous emission, at t = tff and r = r 1 ^, of an 
infinitesimally thin shell moving with the Lorentz factor 
Then one can obtain the formula of the observed flux from the 
jth subjet with viewing angle 9f at frequency v and time T, 



F«>(7> 



2(1 +z)r^cA {i) Acj) iJ, (T)f J) [(l +z)^7 (J) (l - [3 (i > cos 9(T))] 
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(2) 

where z and di are the redshift and the luminosity distance 
of the source, respectively; f U) {v') and A u) represent the spec- 
tral shape and the amplitude of the emission in the comoving 
frame, respectively. Here T = is chosen as the time of arrival 
at the observer of a photon emitted at the origin at t = 0. The 
set of points that emit photons observed at a given time T is 
an arc (or a circle). The functions 9(T) and A<p'\T) represent 
an angular radius and a central angle of the arc, respectively: 



(3) 
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< A0£ b and < 9(T) < A9^-9f), 

cos Afl^-cos e'/> cos $(T) 
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(for otherwise). 



Equation (|3} can be rewritten by 
1 cos 9(T)= ' 
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where t*jL = t^-r^/cft'* is the departure time from the central 
engine of the jth subjet. 

The observed sp ectrum of G RBs is well approximated by 
the Band spectrum (Band et al. 1993). In order to have a spec- 
tral shape similar to the Band spectrum, we adopt the follow- 
ing form of the spectrum in the comoving frame, 

, ., ,(W/i/ "<a»-/3»), 

\^7^ l ''y^' ) (<-^)^ ) exp(^'-<) 



(for z//^ w 



(6) 



where v§'\ a B \ and 0$ are the break frequency and the low- 
and high- energy photon index, respectively. 

As a summary, equations 0, @, and (jfj are the ba- 
sic equations to calculate the observed flux from each sub- 
jet, which depends on the following parameters: Of (which is 
determined by i) u \ ip w , $ obs , and (p° hs through equation Q), 

A6, sub> 7 W > ? dep' r > ' a B . fflt> V o "> A ' ; '' and Z- The wh ° le H g ht 

curve from the GRB jet is produced by the superposition of 
the emissions from the subjets. 



3. Ep-Ei S o CORRELATION FOR A SINGLE SUBJET 

Before examining E p -Ej SO correlation for the multiple sub- 
jet model, it is instructive to calculate E p -Ei SO correlation 
when a single subjet is seen off-axis. Using equations ©, 
(0, Q, and © for Af to t = 1 and a given V , we compute the 
peak energy of the time-integrated spectrum measured in the 
cosmological rest frame, E p , and the "bolometric" isotropic- 
equivalent energy, Ef so , integrating over the 1 - 10 4 keV range 
in the cosmological rest frame. Here the superscript s of 
E s iso means "single". We adopt the following subjet param- 
eters: A6» sub = 0.02 rad, 7 = 300, a B = -1, B = -2.5, and 
yhfQ = 350 keV. In Figure ^ we show E p and Ef so (in units 
of 2.8 x 10 3 7rAr 2 ) for < 6 V < 0.1 rad (the solid line). The 
dashed and dot-dashed lines are E? so QA and E[ m ^ 3 , respec- 
tively. We see that for V > A6> su b, as V increases, both E p and 
Ef so decreases. We focus on small Ef so regime. At first the E p - 
Ei SO correlation approaches E p oc E[ W 0A , but for even larger V , 
EpOcE^ 1 / 3 . This behavior is explained as follows. Firstly the 
spectral peak energy scales as E p oc [l-/3cos(0,, - A6* su b)] _1 
because of the Doppler effect, and for large V , E p oc 6 V ~ . 
Next we compute Ef sg by integrating equation (1) over v and 
T, and study its dependence on 9 V . When 9 V is large but E p is 
in the 1 - 10 4 keV range, the integration over u/(l+z) results 
in a constant depending on the Band spectral parameters and 
another Doppler factor [1 - f3cos9(T)]~ 1 . As for the integra- 
tion with respect to T, we change the variable from T to 0(T), 
and obtain: 

f e ' +Ad ^ A#fl)sin0</0 
X A A'.'.... d-/3cos0)3- U) 

For large 9 V , Acj) ~ A6 sab /9 V , so that Ef so oc 6> v _1 ([l-^cos(6' v - 
A0 sub )]- 2 - [1 - /3cos(0 v + A0 sub )]- 2 )' ~ 9~ 1 [1 - Pcos(0 v - 
A6* su b)] 2 oc E p 2 5 . When 9 V is even so large as E p ~ 1 keV, the 
integration over v/{l+z) results in a factor (1 - f3 cos 9) l+ 13 B , 
so that the same calculation gives us Ef so oc 9 v ~ l [1 — /3cos(0 y — 
A0 sub )]^oc£/for/3 B = -2.5. 

For a single subjet, off-axis events obey E p oc Ej SO 0A for 
small E p regime (but E p > 1 keV). The index of the E p -Ei SO 
correlation, a = 0.4, is obtained irrespective of the intrinsic 
subjet parameters A6* sub , 7, fdep, ro> &b, 0b, and A, as can 
be seen in the above derivation. 

4. Ep-E IS0 CORRELATION IN THE MULTIPLE SUBJET MODEL 

Let us perform Monte Carlo simulations to derive E p -Ei SO 
correlation in the multiple subjet model. For simplicity, we 
generate one GRB jet with opening half-angle A6* tot = 0.3 rad 
and random 5000 lines of sight of the observer with < 
$obs < 0.35 rad according to the probability distribution of 
sinz9 b s d#obs d(p ODS . Then, for each observer, we calculate the 
peak energy of the time-integrated spectrum measured in the 
cosmological rest frame, E p , and the "bolometric" isotropic- 
equivalent energy, E, so , integrating over the 1 — 10 4 keV range 
in the cosmological rest frame. The departure time of each 
subjet t^L is assumed to be homogeneously random between 
t = and t = fdur, where fd ur is the active time of the central 
engine measured in its own frame, and fd ur = 20 s is adopted. 
The central engine is assumed to produce N tot = 350 subjets 
following the angular distribution function 



dN 



1. 



< # u < d c , 



(8) 
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where i3b = A0 tot - A0 su b, and -d c = 0.03 rad. This 
corresponds to the universal structured jet model (see 
iRossi. Lazzati. & Reesl2002llZhang & Meszaros 2002a). The 
angular distribution of the subjets in our simulations is shown 
in Figure |2] The solid circle describes each subjet and 
the dashed circle describes the whole jet. The meaning of 
plus sign will be discussed later. We assume that all the 
subjets have the same values of the following parameters: 

A Cb = °- 02 rad ' %> = 300 < *o = 3 -° x 10 ' 4 cm < a B = -!. 
and 0g = -2.5. The intrinsic spectral parameter ^hv'^ and 
the amplitude A'" are determined so that the time-averaged 
emission from a single subjet viewed on-axis satisfies the fol- 
lowing correlation, 

— — — = e( ESp V (9) 

10 52 erg s" 1 ? V 1 keV / ' 

where L s jso is the time-averaged "bolometric" isotropic - 
equivalent luminosity and E s p is the time-averaged rest- 
frame spectral peak energy of the on-axis emission from 
a single subje t. As for th e validity of this correla- 
tion, [Liang. Dai. & Wu (2004) argue that for long bright 
BATSE GRBs the observed 7-ray flux F is correlated with 
the observed time-resolved £° bs at each time in a simi- 
lar way, i.e., F oc (£'° bs ) 2 , which supports the assumption 
that the on-axis emission of each subjet obeys this corre- 
lation. Llovd-Ronning & Ramirez-Ruiz (2002) show that 
there is a positive correlation between 7-ray luminosity 
and time-resolved rest-frame spectral peak energy by us- 
ing variability-luminosity corr elation (see also Yonetoku et al. 
2004; Ghirlanda et al. 2005b). This correlation could be ob- 
tained by standard synchrotron internal shock model (e.g., 
Zhang & Meszaros 2002b). However, the coefficient £ is 
highly uncertain. Therefore we chose the values of £ so that 
the results of simulations reproduce the observations. We con- 
sider two cases of 7/11/^°' and A 1 ": Case (i) jhi/' Q u ' and £ are 
fixed as 350 keV and 6.0 x 10~ 5 , respectively, for all j. Case 
(ii) ^hv^ and £ are distributed around the above values. 

4.1. Case (i) 

Let us consider the Case (i) as a simple toy model, in 
which all the subjets have the same intrinsic parameters, so 
that we can investigate the pure kinematical effects from the 
multiple discrete emission patches. The results are shown 
in Figure |3 The black solid line shows the E p -Ei SO corre- 
lation for a single subjet derived with the same parameters. 
We see that the black solid line traces the left-side edge of 
the distribution of the simulated bursts. When a single sub- 
jet is seen on-axis, the time-averaged spectral peak energy 
E p = g n hv' a /^{\-(5) « 2g n jhvQ ~ 500 keV, where a numerical 
factor g n (~ 0.7) comes from the contribution of soft emission 
from the whole subjet, while g„ = 1 in the case of point source 
approximation. The observed pulse has a duration determined 
by the angular spreading time as ST = r^AO^ 2 /2c = 2 s. 
Then, according to equation (|9j, L s iso ~ 1.5 x 10 53 ergs s -1 , 
so that E- so = L s iso ST ~ 3 x 10 53 ergs. This corresponds to E p 
reaching its maximum around Ej S0 ~ 3 x 10 53 ergs. When 
more than 1 subjet are seen on-axis, i.e., n s > 2, E p is the 
same as in the case of n s = 1, but Ei S0 ~ n s Ef so . The max- 
imum value of multiplicity n s is about 30, when the line of 
sight is along the center of the whole jet. Then Ei S0 takes the 
maximum value of ~ 10 55 ergs. Points with E p < 500 keV 
correspond to the case of n s = 0, in which all the subjets are 



seen off-axis, i.e., Of > A0 su b for all j. For each line of sight, 
the observed flux is dominated by the emission of the subjets 
with small Of. Thus E p is determined by the minimum value 
of Of, 6>™ n . Let nf f be the number of the subjets with Of 
around 6>™ n . When n° ff = 1, the observed flux is dominated 
by a single subjet, and the #™ n -dependence of E p and Ei S0 is 
determined as discussed in § [3] Such points are on the black 
solid line. When nf > 2, for each 6>™ n , E p is the same as 
for the case of n° ff = 1, but Ei S0 ~ n° s Ef s0 , Thus the scatter 
of the simulated points for E p < 500 keV arises from that of 
n° ff . We find that the right-side edge of the distribution of the 
points follows E p oc £ , 1jo 1//2 . The reason for this behavior is 
as follows. For each 0f n , E p oc [1 - /3cos{0f n - A6> sub )] _1 . 
The other quantity Ej SO is given for the largest n° ff . Since the 
probability that these n° ff subjets have the same axis (fl 01 , ip U) ) 
is quite low, they should be smoothly distributed around the 
line of sight. Then in calculating Ei SO by equation Q for 
the multiple subjets case, we can take A0 ~ tt. Therefore, 
for each 0f n , E iso oc [l-/3cos(6»™ n - A6» sub )]" 2 , and then we 
obtain Ej SO oc E p 2 . Such situation resembles the case of the 
annulus jet model in which the line of sight is inside the 
annul us and the inner ra dius of the annulus changes (see 
Eichler & Levinson 2004). 

For a multiple subjet model, off-axis events (with n s = 0) 
follows E p oc Ei SO a with 0.4 < a < 0.5. This range of a is 
obtained irrespective of the intrinsic subjet parameters A6> su b, 
1, fdep, r , a B , Pb, Vq, and A. 

4.2. Case ( ii) 

We here assume that ^hv'^ is distributed rando mly accord- 
ing to a lognormal distribution function ( Ioka &"Nakamural 
120021) with an average of log(350 keV) and a logarithmic 
variance of 0.2. For given ^hv^, A'" is determined by 
equation (|9j- The coefficient £ is also assumed to obey 
a lognormal distribution with an average of -5 + log(6.0) 
and a logarithmic variance of 0.15. The other parameters 
of the subjets are fixed to the same values as in the pre- 
vious simulation. We calculate Ei S0 and E p , and then as- 
sign a redshift for each observer to calculate the distance 
and the observed lightcurve. The source redshift distribu- 
tion is assumed to be in proportion to t he cosmic star for- 
mation rate. We adopt the model SF2 in Porciani & Madau 
(2001), in which we take the standard cosmological parame- 
ters of H = 70 km s" 1 Mpc" 1 , n M = 0.3, and n A = 0.7. Fi- 
nally, we select detectable events with observed peak photon 
fluxes in the 1 - 10 4 keV band larger than 1.0 ph cm" 2 s" 1 , 
which corresponds to the t hreshold sensitivity of HETE-2 (see 
iBand et aUl2003t lLamb. Donaghv. & Grazianil 12005). Fig- 
ure ^ shows the result of our simulation. Plus signs repre- 
sent bursts that can be detected by HETE-2, while crosses 
represent ones that cannot be detected. They are com- 
pared with the Bepp oSAX and HETE-2 data (points with er - 
ror bars) taken from Ghirlanda, Ghisellini, & Lazzati (2004). 
The solid line represents the best fitted line for 442 GRBs 
with redshifts estimated by the lag-lu minosity correlation 
( Ghirlanda, Ghisellini, & Firmani 2005a). We see that our 
simulated GRBs cover the observed GRBs over three orders 
of E p , so that our multiple subjet model with the intrinsic cor- 
relation E p oc L s im l l 2 under the universal structured jet model 
is consistent with the observations. 

The n s = 1 bursts directly reflect the assumed correlation of 
Ep oc L s ho x l 2 . For larger n s , Ej SO becomes larger, and E p is de- 
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termined by the subjet emission with the largest L s iso observed. 
As a result, all the simulated bursts roughly obeys E p oc Ei SO x l 2 
over about three orders of magnitude in E p . The scatter comes 
from the differences of the number of the observed subjets and 
the differences of the parameters of each subjet. 

5. DISCUSSION 

We have investigated E p -Ei SO correlation in a multiple sub- 
jet model for GRBs, XRRs, and XRFs. We find that off- 
axis events (with n s = 0) for multiple discrete emission re- 
gions show E p oc Ej SO a with 0.4 < a < 0.5. It is assumed that 
the subjet parameters jhvQ and A 01 are distributed so that 
emission of the subjets viewed on-axis follows the correla- 
tion oc U jso 1 , with narrow E s p range (one order of magni- 
tude). Then the Amati correlation (E p oc Ei SO x l 2 ) is reproduced 
over three orders of magnitude in E p . Although the scatter 
around the Amati correlation is large in the simulation, the re- 
sults are consistent with the observed properties of GRBs with 
known redshifts and the BATSE GRBs with pseudo redshifts 
derived from the lag-luminosity correlation. We argue that 
for brighter bursts the Amati correlation arises from intrinsic 
property, while for dimmer bursts it arises from the off-axis 
effects of multiple emissions. The intrinsi c Ej oc L S ,,J^ 2 cor- 
relation is supported by the observations dLiang. Dai. & Wul 
120041: iLlovd-Ronning & Ramirez-Ruizl 120021) and could be 
derived in the context of standard synchrotron internal shock 
model. 

HETE team defines XRRs and XRFs as those events for 
which log[S y (2- 30 keV) /S^(3 0-400 keV)l > -0.5 and 0.0, 
respectively (Lamb et al. 2004). We calculate the observed 
fluence ratio for simulated bursts surviving the peak flux trun- 
cation, and classify them into GRBs, XRRs, and XRFs. The 
ratio of the simulated event rate is Rcrb '■ Rxrr '■ Rxrf ~ 4 : 
3:1. HETE-2 obse rvations show s imilar number of GRBs, 
XRRs, and XRFs (Saka moto et alJl2005l) . We can say that 
the event rate among GRBs, XRRs and XRFs is consistent 
with the observations. Figure [5] shows the redshift distribu- 
tion of GRBs (the solid line), XRRs (the dashed line) and 
XRFs (the dot-dashed line). The mean redshifts of GRBs, 
XRRs and XRFs are 1.9, 3.2, and 2.3, respectively. XRRs 
have a little larger redshifts than GRBs and XRFs. Figure |2] 
plots the viewing angles for detectable XRFs, which are rep- 
resented by plus signs. We see that the main population of 
the XRFs arises from the off-axis effects. On the other hand, 
many XRRs are on-axis events. Since E p > 200 keV for on- 
axis events in our simulation, the on-axis XRRs arise from the 
cosmological redshift effect. The ratio of the on-axis and off- 
axis XRRs is ~ 1 : 1 . We expect that the event rate ratio from 
larger observed samples will give us some information about 
the angular distribution of the subjets within the whole GRB 
jet and the redshift distribution of the GRB sources. 

In this paper, we have performed the simulations with fixed 
Lorentz factor of the subjets, 7 = 300. As discussed in § 
the range of the index a of the E p -Ej S0 correlation for off-axis 
events is independent of the Lorentz factor. We perform the 
same simulations in the case (ii) for 7 = 100 and 500, and 
obtain the Amati correlation (E p oc £ , , so 1//2 ) through all bursts. 
However, the peak photon flux of the XRF is small for lower 
Lorentz factor. For 7 = 100, we obtain Rcrb ■ Rxrr ■ Rxrf ~ 
15 : 10 : 1. Alternatively for 7 = 500, R C rb ■ Rxrr ■ Rxrf ~ 3 
: 2: 1. 

Figure [6] shows the distribution of the durations in the 
50-300 keV band for GRBs (solid line), XRRs (dashed line), 



and XRFs (dot-dashed line). GRBs have a bimodal distri- 
bution as observed by BATSE. We have already shown why 
GRBs have the bi modal duration d istributi on in our multi- 
ple subjet model (Toma, Yamazaki, & Nakamura 2005): the 
Tgo duration of an n s = 1 burst is determined by the width 
of a single pulse, while that of an n s > 2 burst is deter- 
mined by the time interval between the observed first pulse 
and the last one. These two different timescales naturally 
lead to a division of the burst Tgo durations into the short 
and long ones. We also calculate the distribution of the 790 
durations in the 2-25 keV band. Figure [7] shows the re- 
sult. These distributi ons are not inconsistent with the HETE- 
2 data (see Fig. 4 of Sakam oto et alJ l2005). The Tgo dura- 
tions of n s = 1 bursts (i.e., short bursts) become larger when 
they are measured in the softer band, since soft emission 
from the perip hery of the subjet is observ ed for a longer 
time. Yamazaki, Ioka, & Nakamura (2004b) have predicted 
short XRRs in our unified model, which are confirmed in 
this simulation. These are events of a single subjet viewed 
off-axis or viewed on-axis with slightly high redshift. In- 
deed, GRB 040924 may be an example of short XRRs, from 
whic h recent HST o bservation reveals the supernova signa- 
ture (ISoderberg et al.ll2.005h . This event supports our unified 
picture. 

In this paper, we considered the #~ 2 -angular distribu- 
tion of the subjets. Averaging by a solid angle satisfy- 
ing (A# su b) 2 < VL < (A# t ot) 2 , the distribution of the emis- 
sion energy (or almost equivalently the angle-averaged ki- 
netic energy) is the same as the universal structured j et model 
jRossi. Lazzati. & Reesl 120021 Sang & M eszarosl l2002al) . 
The universal structured jet model has been criticized by 
Lamb, Donaghv, & Graziani (2005): in the universal struc- 
tured jet model, it is assumed that XRFs are observed when 
the jet is viewed from fairly large angle, so that the model 
overpredicts the number of XRFs, which is inconsistent with 
the observed ratio of the number of XRFs a nd GRBs de- 
tected by HETE-2. Then, IZhang et al.1 J2004I) modified the 
universal structured jet model, and showed that if the jet 
is structured with a Gaussian-like shape, the number of 
XRFs becomes small. In these works it is assumed that 
the jet is continuous and th ere are no cold spots insi de the 
jet. As shown in this pap er, Eichl er & Levinsonl (120041) . and 
Yamazaki. Ioka. & Nakamura (2004b), if the observer points 
toward the cold spot (i.e., n s = 0), XRFs or XRRs are ob- 
served. While if n s > 2, the event looks like a long GRB 
irrespective of the viewing angle. In our model, the ratio 
of the total solid angle with n s > 2 and n s = determines 
the event rate of GRBs and XRRs/XRFs. Interestingly, we 
find that the power-law profile with an index of -2 is prefer- 
able to the Gaussian profile in order to reproduce the ratio 
of observed event rate of GRBs, XRRs, and XRFs, because 
the solid angle with n s = is small in the Gaussian profile. 
lLazzati & Begelmanl (120051) have recently argued that in the 
context of the collapsar model, 6~ 2 angular profile might be 
obtained as a consequence of the physics in the jet breakout 
irrespective of the jet structure inside the progenitor. From the 
observational side, we can estimate the pseudo jet opening an- 
gle distribution. Using the Ghirlanda correlation (E p oc £ 7 ' 71 ) 
where E 1 = Ei SO B 2 /2 (Ghirlanda, Ghisellini, & Lazzati 2004) 
and the Yo netoku correlation (E n oc £„ 05 ) where L„ is ap eak 
luminosity JYonetoku et al.l2004l) . lYonetoku et"aTl (120051) ob- 
tained that the pseudo jet opening angle obeys f(0j)d0j oc 
9~ 2 d9j. This is compatible with the power-law structured jet 
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model: if all bursts were observable, the distribution would be 
uniform per unit solid angle and f(6) oc 6. However £„„ for 
the smaller viewing angle is brighter by a factor of 9~ 2 , so that 
the maximum observable distance is larger by a factor of 9~ l 
which contains a volume larger by a factor of 6~ 3 . Then we 
have/(6»)ocr 2 . 

Late phase evolution of a set of multiple subjets is rather 
complicated and hard to be predicted. Cold spots do not 
produce high energy emission but may be filled with the ki- 
netic energy that is not d issipated at small radius (see also 
iLevinson & Eichlerl i2005). Even if cold spots are not filled 
with the kinetic energy, all subjets begin to expand sideways 
and would merge into one shell. In any case, late afterglow be- 
havior may be well approximated b y the results from the con- 
tinuous structured jet model (e.g., Kumar & Granot 2004J). As 
shown in Fig. [2 almost all XRFs arise when all the subjets are 
viewed off-axis, i.e., n s = 0, while the observers see the whole 
jet on-axis. Then, the late phase (> 1 day) properties of XRF 



Amati, L. et al. 2002, A&A, 390, 81 
Amati, L. et al. 2004, A&A, 426, 415 
Band, D. L., et al. 1993, ApJ, 413, 281 
Band, D. L. 2003, ApJ, 588, 945 
Butler, N. R., et al. 2005, ApJ, 621, 884 
Eichler, D., & Levinson, A. 2004, ApJ, 614, L13 
Fynbo, J. P. U., et al. 2004, ApJ, 609, 962 
Ghirlanda, G., Ghisellini, G, & Lazzati, D. 2004, ApJ, 616, 331 
Ghirlanda, G, Ghisellini, G, & Firmani, C. 2005a, MNRAS, 361, L10 
Ghirlanda, G, Ghisellini, G, Firmani, C, Celloti, A., & Bosnjak, Z. 2005b, 
MNRAS, 360, L45 

Granot, J., Ramirez-Ruiz, E., & Perna, R., 2005, preprint I astro-ph/0502300 i 

Ioka, K., & Nakamura T. 2001, ApJ, 554, L163 

Ioka, K., & Nakamura, T. 2002, ApJ, 570, L21 

Kumar, P., & Piran, T. 2000, ApJ, 535, 152 

Kumar, P. & Granot, J. 2004, ApJ, 591, 1075 

Lamb, D. Q., et al. 2004, NewA Rev., 48, 423 

Lamb, D. Q., Donaghy, T. Q., & Graziani, C. 2005, ApJ, 620, 355 

Lazzati, D. & Begelman, M. C. 2005, preprint i astro-ph/0502084 i 

Levinson, A. & Eichler, D. 2005, preprint i astro-ph/0504125 i 

Liang, E. W., Dai Z. G, & Wu, X. F. 2004, ApJ, 606, L29 

Lloyd-Ronning, N. M., & Ramirez-Ruiz, E. 2002, ApJ, 576, 101 



afterglows may not be like orphan afterglows but may show 
simila r behavior to those of normal GRBs (e.g., Am ati et alJ 
2004). On the other hand, as rare cases, when the whole 
jet is viewed off-axis, XRF aftergl ows may resemble the or - 
phan afterglow (e.g.. iGranot. Ramirez-R uiz. & Perna 2005). 
XRF 030723 may be a member of such a class iButler et alJ 
l200llFvnbo et al.l2004l) . 
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Angle [rad] 

FIG. 2. — Angular distribution of Mot = 350 subjets confined in the whole GRB jet in our simulation. Each subjet is located according to the power-law 
distribution function of eq.j8}. The whole jet has an opening half-angle of A£> to t = 0.3 rad. The subjets have the same opening half-angles of A6> SUD = 0.02 rad. 
The angular size of the subjets are represented by the solid circles, while the whole jet is represented by the dashed circle. The viewing angles for detectable 
XRFs in our simulation are represented by plus signs. 
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FIG. 4. — Same as Figure[3] but in the multiple subjet model in which subjet parameters 7^q'*' ) , A {i \ and § are distributed (see text for details). Plus signs 
represent bursts that can be detected by HETE-2, while crosses represent o nes that cannot be detected. They are compared with the BeppoSAX and HETE-2 data 
(points with error bars) taken from Ghirlanda. Ghisellini, & Lazzati 1 2004). The solid line represents the best fitted line for 442 GRBs with redshifts estimated 
by the lag-luminosity correlation i Ghirlanda. Ghisellini, & Firmani 2005a). 
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FIG. 5. — Redshift distribution of the simulated bursts surviving the peak photon flux truncation. Solid line, dashed line, and dot-dashed line represents the 
distribution for GRBs, XRRs, and XRFs, respectively. 
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FIG. 6. — Distribution of the T 90 durations in the 50-300 keV band of the simulated bursts surviving the peak flux truncation. Solid line, dashed line, and 
dot-dashed line represents the distribution for GRBs, XRRs, and XRFs, respectively. 




FIG. 7. — Distribution of the T 90 durations in the 2-25 keV band of the simulated bursts surviving the peak flux truncation. Solid line, dashed line, and 
dot-dashed line represents the distribution for GRBs, XRRs, and XRFs, respectively. 



